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CANDIDATES OF HALPHA EMITTING REGIONS IN MAGELLANIC STREAM IV CLOUD 
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ABSTRACT 

From Ha narrow band observations, we identified three Ha emitting regions in the direction of Magellanic 
Stream IV (MS IV). They consist of three parallel filaments with 2 arcmin width and 6-30 arcmin length at 
12 arcmin intervals. The mean surface brightness of them is ~ 2 x 10~^^ erg s~^ cm~^ arcsec"^. Because 
of their low surface brightness, the regions were not detected in previous Ha surveys. In H I map, the position 
of the filaments overlap MS, suggesting that they are parts of MS, but there also exists a local Hi structure. If 
the filaments associate with MS, the sizes are 30 pc x 100 - 500 pc. The filaments lie at the leading edge of a 
downstream cloud, which supports a shock heating and its propagation (shock cascade) model for the ionizing 
source. If they are local objects, on the other hand. Fossil Stromgren Trails of more than two stars is a possible 
interpretation, and the sizes would be 0.1 pc x 0.3 - 1.5 pc at 180 pc distance. The positional information of the 
Ha filaments presented in this letter enables us future spectroscopic observations to clarify their nature. 

Subject headings: galaxies: individual (Magellanic Stream) — intergalactic medium — Galaxy: halo 



1. INTRODUCTION 

Magellanic Stream (MS hereafter) was found as a long fil - 
ament of H I gas around our Galaxy ( [Mathewson et al.||1974| ). 
MS is thought to be stripped gas from Small Magellanic Cloud 
and/or Large Magellanic Cloud. It has been revealed that MS 
consists of multi-phase gas component s; molecular g as, H I gas, 
H II gas, and highly ionized plasma (F ox et al.|2010 and refer- 
ences therein). MS is the only clear example of a gaseous halo 
stream in the close proximity of Milky Way ( [Stanimirorovic et 



al. 2008). It gives us important clues to the physical interaction 
between Galaxy and its satellites in detail. 

The Ha emission from MS has been first reported bylWeiner 



In order to study the ionization mechanism and kinematics of 
the MS Ha clouds, detailed spectroscopic follow-up observa- 
tions are required. However, positional information necessary 
to perform a medium to high dispersion spectroscopy has been 
so far very limited. The exact position and detailed morphol- 
ogy of the Ha emitting region found by WW96 was not clear, 
because they used a Fabry-Perot interferometer and the spatial 
structure inside the Fabry-Perot beam of 7 arcmin diameter was 
averaged. 

There are several other Ha surveys which overl aps the MS 
fields. S outhern Ha Sky Survey Atla j^ (SHASS A; ' Gaustad'et 



|& WilHamsl ( [T9961 hereafter WW96). |Putman et al.] ( |2003b| ) 
have found several other Ha emitting regions. The ionizing 
source of these regions is still a mystery because no ionizing 
star is observed in the stream. There are several models to ex- 
plain the origin of the MS Ha clouds. WW96 argued that the 
emissions are explained by ram pressure heating. However, this 
model required much higher density of halo gas than WW96 es- 
timated, because the distance to MS was found to be larger than 
expected ( [Weiner et al. 2001 ). Bland-Hawthorn & Maloney| 
( |1999| ) presented another model that ionizing photons escaping 
from Galactic disk heat up MS. This model predicted that the 
Ha emission would be smoothly distributed. However it is in- 
consistent with the fact t hat the contrast of t h e Ha emission is 
as large as a factor of 30 ( [Weiner et al.|200ll ). |Bland-Hawthorn 



|al.|2001| ) presents Ha map of 0.8 arcmin pixel. They used a 3.2 
nmjwidth filter to detect Ha emissions. Wisco nsin Ha map 



|et al.| ( |2007 , hereafter B07) presented yet another model incor- 
porating a shock heating and its propagation (shock cascade). 
They predicted high Ha/H/3 ratio, and strong [O ll], [S ll] emis- 



sions, which have not been examined yet. Recently, |Fox et 
al. ( 2010| ) investigated the absorption features of active galac 



tic nuclei behind MS in order to study the ionization structure 
of the MS Ha clouds. They found that the ionization struc- 
ture cannot be reproduced by an evaporative interaction with 
the Galactic corona plus photoionization. 



pe]0(WHAM; Raynolds, et al.|l998[|Haffi^er^ al.,2003 ) pro- 
vides another Ha map with high velocity resolution in a veloc- 
ity range from —80 km s~^ to +80 km s~^ with respect to the 
local standard of rest (LSR) with 1 degree beam using a Fabry- 
Perot spectrometer. The spatial resolutions of those works are 
still not enough for a slit spectroscopy. 

We performed imaging observations with higher spatial res- 
olution at the northern part of MS IV where WW96 detected 
the Ha emissions. In this letter, we report the positions and the 
shapes of the Ha emitting regions we detected. 

2. OBSERVATION 

We observed two WW96 fields, MS IV-C (23:44:34.8,- 
12:23:20,12000.0) an d MS IV-D (23:4 0:49.0,-10:45:32,12000.0) 
with 2kCCD camera ( [Itohet al.|2001| ) mounted on Kiso 105 cm 
Schmidt camera at the Kiso Observatory in September 2011. 
The pixel scale is 1.5 arcsec pixel" ^, and the field of view 
(FoV) is 51.2 arcmin square. We used R-band and three narrow 
band filters (N657, Ha6577, and Ha6417). The filter charac- 
teristics are shown in Table 1 . The variation of the wavelength 
coverage is no worse than 3.6A across the filter in Ha6577 
and Ha6417. For N657, the transmission curve at the corner is 
shifted by lOA to the red from that at the center. We took the 
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exposures with ~ 7.5 arcmin dithering. The typical seeing size 
was 4.1 arcsec. 

The data reduction was performed using the mosaic CCD 
data reduction packages (Yag i et al._2QQ2 K Overscan subtrac- 
tion, bias pattern subtraction, flat-fielding with domeflat were 
performed in a standard manner. The pixel scale of the im- 
ages are slightly different among filters. The difference was 
corrected before coadding using the results by SMOKA team 
(S. Ichikawa, private communication) [j The sky background 
was subtracted in each frame before coadding with background 
mesh size of 1024 pixels (25.6 arcmin). Therefore, structures 
larger than 25 arcmin square would have been removed. Af- 
ter the coadding, stellar objects were detected by SExtractor 
( [Bertin & Amouts||l996| , and the astr ometric calibration was 
performed using wcstools ( Mink|2002| ) referring to USNO-Bl 
catalog ( |Monet et al. r2003'). 

We calibrated the flux using the photometry of stars by Sloan 
digital sky survey (SDSS) data release 7 ( [Abazajian et al.|2009 ) 
in the northern half of the MS IV-D field. The color conversion 
from SDSS r and i to Kiso 2kCCD filters are constructed by 
fitting a quadratic function to the photometry of model spectral 
energy distributions (SEDs) by Bruzual-Pearson-Gunn-Stryker 
atlas [^ and the filter responses, as the method in the appendix of 
Yagi et al. ( 2Q10| ). The best-fit functions are 



Kiso ■■ 



' + Co + ci(r - i) + C2(r - i^ 



(1) 



where (co,ci,C2)= (0.003,0.267,0.052), (-0.075,0.773,-0.292), 
(-0.034,0.653,-0.210), and (0.014,0.262,0) for R, N657, 
Ha6577, and Ha6417, respectively, in AB -magnitude for 
—0.1 < r — i < 1.0 stars. The dispersion of our photome- 
try around the fit is 0.06 - 0.1 mag in rms converted from the 
median of absolute deviation in each band. As MS IV-C is not 
covered by SDSS, we roughly estimated the zero point using R- 
band (Rp-band) photometry of USNO-Bl. ( |Monetetal.|2003|) 
in Naval Observatory Merged Astrometric Dataset ( |Zacharias| 
et al.|2 005). The error is estimated to be ~ 0.15 mag. Galac- 
tic extincti on is about Ar = 0.0 8 - 0.09 in MS IV-C and MS 
IV-D from Schlegel et al. ( 1998| ) via the extinction calculator in 
NASA/IR\C extragalactic database (NEdQ We adopted 0.09 
mag for all of the 4 bands. The limiting surface brightness (SB) 
of the final images of MS IV-D in 15 arcsec (10 pix) square are 
presented in Table 1 . 

3. DETECTION OF THREE Ha EMITTING REGIONS 

We performed eye inspections to detect possible Ha emitting 
regions. In the MS IV-C field, we could not find any sign of Ha 
excess brighter than 800 mRJjwith size between 15 arcsec and 
25 arcmin scale in the 54 x 54 arcmin square field. The conver- 
sion from SB of N657-R to total Ha+[Nll] brightness requires 
the redshift and the line full width at half maximum (FWHM) 
information, which are unavailable yet. In this study, we there- 
fore assume that the line is located at the wavelength of the 
maximum transmittance of each filter, and that the line width 
is much narrower than the filter FWHM. This gives the lower 
limit of the Ha+[Nll] brightness. 

In the MS IV-D field, we detected regions with N657 excess 
relative to R-band. Fig. la and lb show the R-band subtracted 
N657 image (N657-R image) and the N657 image of the MS 
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IV-D field, respectively. The images were binned by 10x10 
pixels (15 arcsec square) to increase signal to noise ratio (S/N). 
Three filamentary Ha regions are seen in the lower half region 
of Fig. la and lb. We call the Ha regions as Dl, D2, and D3 
from the east to the west, hereafter. Except the three, we could 
not detect Ha emitting structures brighter than 240 mR (la of 
the N657-R image) with size between 15 arcsec and 25 arcmin 
scale in 74 x 60 arcmin field of MV IV-D. Fig. Ic shows the 
effective exposure time map of N657-R and the position of the 
Ha regions. Because of the bad weather, the southwest part has 
low S/N. It should also be noted that a bright star at the south- 
west of D2 made it difficult to measure the emission around 
it. 

We found Dl and part of D3 in the Ha6577 image with SB 
of 26.3 mag arcsec"^. The three Ha regions are fainter than 
27 mag arcsec"^ in R-band. In the Ha6417 image, neither Dl 
nor D2 were detected, and D3 was located outside of the FoV 
Dl was also found in the N657-Ha6417 image, and marginally 
recognized in the Ha6577-Ha6417 image. In the Ha6577-R im- 
age, Dl and D3 were detected, and D2 was marginally detected. 
In the N657-Ha6577 image, Dl was marginally positive (N657 
bright) with low S/N, and D2 and D3 were not detected. In 
summary, the excess was found in N657-R, N657-Ha6417, and 
Ha6577-R, and marginally in N657-Ha6577. The detection in 
several combination of bands suggests that the features are real 
emissions. 

We investigated the possibility of artifacts. If the filaments 
were fringes made by sky lines, they should also appear in other 
fields, but we could not detect such features in other fields (e.g., 
MS IV-C). We also checked whether the regions were seen in 
coadded images of each day, and confirmed that they were not 
moving objects nor transient phenomena. The positions and 
the shapes remained the same during the observation. When 
we observe other brighter stars in another field, we found some 
ghosts of the optics appear. The ghosts moved when we shifted 
the FoV, and they appeared at different position in Ha6577 and 
N657. On the contrary, D1-D3 did not move and stayed at the 
same position in the sky in both filters. Moreover, we made ad- 
ditional observations in November 2011 for the filaments with 
Ha6577. At that time, we shifted the field to avoid the bright 
star at south-east. The filaments Dl and D3 were still detected 
at the same position. D1-D3 are therefore not ghosts. We con- 
clude that the three regions are genuine Ha emitting regions. 

The flux weighted center of the regions are shown in Table 
2. The positions are 20-50 arcmin away from the aperture of 
MV IV-D by WW96 (Fig. Ic). 

We defined the boundary of the filaments for a quantitative 
analysis as follows. First, we masked bright objects in R-band, 
because a residual of subtraction may remain in the N657-R 
image at the position. Next, we applied a 7x7 pix (10.5 arcsec 
square) median filtering to gain S/N. Then the R-band image 
is subtracted from the N657 image, and binned into lOx 10 pix 
(15 arcsec square). In the image, we defined the boundary as 
the isophote of 220mR, which corresponds to l-a. Above the 
isophote the regions have ~2 arcmin wide and 16, 6, and 29 
arcmin long for Dl, D2, and D3, respectively. The interval be- 
tween Dl and D2 are about 12 arcmin, and that between D2 
and D3 is also ~12 arcmin. If the regions lie at MS distance 



' erg s ^ cm ^ arcsec ^ at 6563 A, and 1 R=1000 mR. 
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(55 kpc; Putman et al.||20Q3? ), the width of the three Ha re- 
gions is 30 pc, intervals between them are 200 pc each, and the 
length are 100- 500 pc. 

For photometric analysis, we used the isophote as the bound- 
ary and extracted the regions before median filtering. The SB of 
the Ha regions in the N657-R image are 370 mR on average in 
15 arcsec square, which corresponds to Ha+[N ll] brightness of 
~ 2.1 X 10~^^ erg s~^ cm~^ arcsec"^. The SB gradient is seen 
in Dl and D3. The north-west edge is the faintest in SB in D3, 
and the gradient is about 20 mR arcmin"^. The gradient in Dl 
has an opposite direction; the SB is the faintest at the southeast 
edge, rises by 60 mR arcmin"^ for 10 arcmin, then gradually 
get fainter toward north west. The gradient is not seen in D2. 

4. COMPARISON WITH Ha SURVEYS 

Our observe d field was cov ered by two Ha surveys. One is 
SHASSA ( Gaustad et al.|[2Q0 1 ). The H a bandwidth (3.2n m) 
of SHASSA is comparable to our N657. |Finkbeiner| ( [25Q3] ) re- 
moved stellar contamination from the SHASSA data and cal- 
ibrated the flu^l^ There are several bright parts in SHASSA 
data which overlap D1-D3. However, when we checked the 
original data of SHASSA, they are not permanent but variable 
among exposures. The bright parts would be artifacts. 

The other survey is WHAM, which covers a velocity range 
of ^ ±80km s"^ with respect to LSR. One of WHAM be am, 
(/, b) =(13311, -67.054) is just on Dl ( [Haffner et al!pQQ3) . If 
we average the Ha+[N ll] flux of the regions D1-D3 in 1 degree 
aperture of WHAM, it would be observed as '^ 5 mR excess, 
which is too faint to be found in WHAM. 

5. COMPARISON WITH H I MAP 

We overlap the position of the Ha regions on the H I contour 
using Galactic All-Sky Survey (GASS) Second Data Release 



( McClu re-Griffiths et al.||20Q9| [Kalberla et aLl |2010p1in Fig. 
2. The regions D1-D3 lie on an overlap of two H I components 
at vlsr ^ —210 km s~^ and vlsr ^ —30 km s~^. We thus 

coadd the GASS data over two velocity ranges; —170 250 

km s~^ as the high velocity component(HVC), and —15 63 

km s~^ as the low velocity component(LVC). 

In HVC map (Fig 2a), an Hi cloud wit h a horseshoe lik e 
morphology is located at (/, b) =(73.5, -68) ( [Brtins et al.|2005| ), 
which is at the southeast of our observed field. The H I clo ud is 
also recognized in the —211.2 km s~^ panel of Fig. 6 of Put- 



manetal.| ( |2003a| ). They mentioned that filaments of MS appear 
to merge at several points, and there are dense concentrations 
of gas at these points. They listed the cloud at (l,b)=(74,-68) 
as an example of the merging points. The regions D1-D3 are 
located on this cloud. Spatial correspondence between D1-D3 
and the H I cloud suggests that the Ha regions are associated 
with MS, and we regard them as candidates of Ha emitting re- 
gions of MS. The H I column density around the Ha regions are 
2 - 6 X 10^^ cm"^ from the contour by [Brtins et al.| ( |2005| ). 

In LVC (Fig. 2b), there are two clouds around the edge of 
the filaments, and it is also possible that D1-D3 are associated 
with LVC. The H I column density around the Ha regions are 
3-6 X lO^^cm-^. 

6. COLUMN DENSITY OF THE IONIZED GAS 

We roughly estimate the column density of the ionized gas. 
As the typical flux is 370 mR including [Nil], we adopt 300 

^^ http://skymaps.info/ 

^ ^ http://www.astro.uni-bonn.de/hisurvey/ 



mR for Ha flux in this section, assuming [Nll]/Ha ~ 0.2 ( [l^ 
et al.||2007| ) for simplicity. The value is a typical one for H II 
regions, and it would be largely different if the ionization pro- 
cess is different, such as shocks. We assume that the Ha emit- 
ting region is a homogeneously distributed, optically thin gas 
cloud, and is embedded in an isotropic radiation field. For Ha, 

1 Ray lei^h corresponds to an emission me asure (EM) of about 

2 cm"^ pc( |McCunough & Benjam in'2001 ). Therefore the EM 
of D1-D3 is ^0.6 cm"'^ pc. The EM is written as 

EM = {nl)L = 0.6 cm"^pc, (2) 

where ng is the electron density, and L is the mean length of 
line of sight. The column density is then 

1/2 

(3) 



^(^L = 7.6 X 10 



18 



lOpc 



cm 



It is about 1/2-1/8 of H I column density (2 - 6 x 10^^ cm-^), 
if L = 10 pc. 

This rough calculation may overestimate the column density, 
because, in the real, the Ha surface brightness distribution is 
not homogeneous. In addition, the redshift and the degree of 
[N II] contamination are unknown. 

7. COMPARISON WITH MODEL BY B07 

Since WW96 discovered the Ha emission from MS, several 
hypotheses have been proposed for the ionization mechanism of 
the Ha emitting gas. A promising theory would be the shock- 
induced model presented by B07. They suggested that the inter- 
action between the fragmented clouds formed in the upstream 
of MS and the tenuous gas ablated from the upstream clouds 
would lead to shock ionization at the leading edges of the down- 
stream clouds. 

If the northwest cloud in Fig. 2a, which lies in the down- 
stream of the main MS IV cloud at southeast, is following the 
main cloud. Ha would be emitted at the leading edge of the 
northwest cloud according to the B07 model. On the contrary, 
D1-D3 appear to lie along the H I stream between the northwest 
and the main clouds. Their straight and almost parallel filamen- 
tary morphology do not resemble the Ha distribution predicted 
by B07. However, note that the widths of D1-D3, ^ 30 pc, 
are smaller than the spatial resolution of the simulations (^40 
pc) of B07. In addition, kpc scale structures of the Ha emis- 
sion may be lost in this study due to our sky subtraction method 
(see Section 2). Thus, the spatial resolution mismatch between 
our observation and the B07 simulation may lead this apparent 
inconsistency in ionized gas morphology. The Ha regions we 
detected may be a part of dense ridges, which can be formed by 
Kelvin-Helmholtz instability at the H I cloud surface, of tenu- 
ously extended ionized clouds. 

8. POSSIBILITY OF FOSSIL STROMGREN TRAIL 

The long and str aight shape of D1-D3 resembles t he Ha fil- 
ament reported by I McCullough & Benjamin] ( |2001[ hereafter 
MBOl). The MBOl cloud has a region of 20 arcsec wide and 
1.2 degree long and forks into two segments with the widest 
separation of 5 arcmin. 

MBOl presented four possible origins for such a cloud; 1) a 
jet, 2) a filamentary local nebula, 3) a shock ionized trail, and 4) 
a Fossil Stromgren Trail (EST), and they concluded their object 
would be EST. If D1-D3 are local objects, the origin of D1-D3 
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would be one of the four. Following the discussion by MBOl, 
the faint SB rejects 1), and the straight shape and the intensity 
gradient along the shape reject 2), and 2 arcmin width rejects 
3) as the D1-D3 origin. And 4), FST by a single star, is dif- 
ficult, since the gradient has opposite direction in Dl and D3. 
If D1-D3 are local objects, the possible origin is that they are 
FSTs of different stars. All the three might be made by three 
different stars or by two stars if D2 is a branch of Dl or D3 like 
the Y-shaped filament found by MBOl. 

If D1-D3 are FSTs, we can estimate the distance using equa- 
tion (5) of MBOl. Putting I(Ha)obs=300 mR, width 7=2 ar- 
cmin, and the assumption by MBOl that projection factor p=l, 
temperature T=10^K, and density n=l cm~^, we can estimate 
that the distance is '^180 pc. The physical size is O.lpc wide 
and 0.3 - 1.5 pc long. 

9. SUMMARY 

We detected three filamentary Ha emitting regions in the MS 
IV-D field. Their sizes are 2 arcmin wide and 6-30 arcmin 



long, and their surface brightness of Ha+[Nll] are 370 mR on 
average in 15 arcsec bin. They are candidates of Ha emitting 
regions associating with MS IV, although we cannot deny the 
possibility that they are local objects. If they are at MS distance, 
the size is 30 pc x 100 - 500 pc. If they are local objects. Fos- 
sil Stromgren Trails ( McCullough & BenjaminpOO l*) created 
by more than two stars is a possible interpretation. Future spec- 
troscopic follow-up will reveal whether they are a part of MS. 
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Fig. 1 . — Coadded images of the MS IV-D field: N657-R (a) and N657 (b). North is up, and east to the left, and the coordinates are J2000.0. The image size is 
87.5 X 100 arcmin. The scale of colorbar is in units of Rayleighs. Effective exposure time map of N657-R is shown as (c). The colorbar is in units of hours. The 
aperture of WW96 (cyan) and the position of the Ha emitting regions (green) are overlaid. 




Fig. 2. — Overlay of N657-R map in this studv with D1-D3 mark on the H I contour. The size is 5 x 4 degree, north is up and east to the left. The H I contour 
is created from GASS Second Data Release (McClure-Griffiths et al. 2009 Kalberla et al. 2()10). (a) The velocity of -170 < v < -250 km s~^ are coadded. 
The contour corresponds to 2x lO-*^^ cm~^ interval from 2x lO^'^ cm""^. (b) The velocity of —63 < v < —15 km s"-*^ are coadded. The contour corresponds to 
1x10^9 cm-2 interval from 3x10^9 cm-^. 
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Table 1 
Used Filters, and 1-a surface brightness(SB) of the MS IV-D field. The column v(Ha) shows the recession velocity 

COVERAGE FOR Ha. THE 1-cr SB IS ESTIMATED AROUND THE HIGHEST FLUX PART IN FiG. IC, AND MEASURED IN 15 ARCSEC SQUARE. 

For N657 and Ha6577, the SB is also given in units of milli-Rayleigh(mR). 



name 


R.A.(J2000) 


Dec.(J2000) 


I 


b 


MS IV-Dl 
MS IV-D2 
MS IV-D3 


23:40:43 
23:40:13 
23:39:17 


-11:09:06 
-11:19:40 
-11:21:00 


73.3891 

72.8502 
72.3947 


-66.9449 
-66.9864 
-66.8430 



Table 2 
The flux weighted center of the Ha emitting regions 



